Fano effect in a Josephson junction with a quantum dot 
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We theoretically investigate the Fano effect in dc Josephson current at the absolute zero of tem- 
perature. The system under consideration is a double-path Josephson junction in which one path is 
through an insulating barrier and the other one is through a quantum dot (QD). Here the Kondo 
temperature is assumed to be much smaller than the superconducting gap, and the Coulomb inter- 
action inside the QD is treated by the Hartree-Fock approximation. It is shown that the Josephson 
critical current exhibits an asymmetric resonance against the QD energy level. This behavior is 
caused by the interference between the two tunneling processes between the superconductors; the 
direct tunneling across the insulating barrier and the resonant one through the QD. Moreover, we 
find that the Josephson critical current changes its sign around the resonance when the Coulomb in- 
teraction is sufficiently strong. Our results suggest that 0-n transition is induced by the cooperation 
of the Fano effect and the Coulomb interaction inside the QD. 

PACS numbers: 73.63.Kv, 74.50.+r. 



I. INTRODUCTION 

In the last decade, transport properties of quantum 
dot (QD) systems have been extensively studied. These 
studies are motivated not only by their potential appli- 
cationsii^ but also by fundamental interests in many 
body problems and interference phenomena, such as the 
Kondo effect 4 .^ and the Fano effec^^aLi 2 . 

Kobayashi et al. observed the Fano effect using a 
QD embedded in an Aharonov-Bohm interferometer—. 
In this system, there are two transmission processes be- 
tween the source and drain electrodes. One is through 
the continuum in the arm, and the other is through the 
discrete energy level of the QD. The interference be- 
tween these processes results in an asymmetric resonance 
in the differential conductance. Thereafter there have 
been many studies on the Fano effect in QD systems, 
but most of these studies have focused on normal metal 
systems^ ' 10 ! 11 . Thus, the Fano effect in superconducting 
system o 13 ! 14 is less understood. 

Meanwhile, many authors have investigated supercon- 
ductor/quantum dot/superconductor (S/QD/S) junc- 
tion^^!^1122^L2^.. In these systems, the ratio 
of the Kondo temperature Tk to the superconducting 
gap A is a key parameter. In the strong coupling limit 
Tk > A, the Kondo effect survives even in the pres- 
ence of the superconductivity; a Cooper pair is broken 
in order to screen the localized spin in the QD. On 
the other hand, in the weak coupling limit Tk <C A, 
the Kondo effect is negligible because a strongly bound 
Cooper pair cannot be broken. Then, the Cooper pair 
feels the localized magnetic moment in the QD. Under 
this situation, when the Coulomb interaction is strong 
inside the QD, so-called 0-n transition occur o 24 i 25 i 26 i 27 i 28 . 
It means that the dependence of the Josephson current 



on the phase difference 9 changes from I = \I C \ sin 9 to 
I = \I C \ s'm(9 + 7r) = — 1 7c | sin#, i.e., the critical current 
becomes negative. 

Recently, Zhang studied the Fano effect in a Josephson 
junction with a QD coupled in parallel 13 . The calcula- 
tions were done within the finite U slave boson mean field 
theory for Tk > A. It was concluded that the Fano effect 
caused 0-ir transition in that regime. 

In this paper, we calculate the Josephson current 
through the similar system as in Ref.dH, but for Tk -C A. 
The purpose of this study is to examine whether the 0- 
7T transition caused by the Fano effect occurs or not in 
this regime. The system we consider is a double-path 
Josephson junction; a S/QD/S junction and a conven- 
tional Josephson junction are connected in parallel (see 
Fig. [T]). Here we employ the Hartree-Fock approxima- 
tion (HFA) in treating the Coulomb interaction inside the 
QD. It is found that the critical current exhibits the char- 
acteristic Fano-like dependence on the QD energy level. 
We also show that the critical current near the resonance 
can change its sign under the strong Coulomb interac- 
tion. This means that the combination of the Fano effect 
and the Coulomb interaction effect causes 0-7r transition. 
Finally, by using a simple perturbative approach with 
respect to the electron tunneling, we discuss the physi- 
cal origin of the transition and estimate the parameter 
regime where the system behaves as a 7r-junction. 

The organization of this paper is as follows. In the 
following section, we introduce the model Hamiltonian 
and present the details of our calculation. In Sec. HI, 
we show the results for the non-interacting case and the 
interacting case respectively. Discussions are given in 
Sec. IV with the results of the perturbation expansion. 
Section IV is devoted to a brief summary. We also present 
the results obtained by the Hubbard-I approximation 29 
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FIG. 1: Schematic illustration of the system under consid- 
eration. A S/QD/S junction and a conventional Josephson 
junction are connected in parallel. 
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Here we introduced the lesser Green functions, which are 
defined by 
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in the Appendix. 



II. MODEL AND FORMALISM 

In order to describe the system shown in Fig. [TJ we 
use the following Hamiltonian: 



W /3 = X! & C 0k,<? C Pk,<r 



(1) 
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W D = ^e d <4<f CT + Unfn L , (3) 

A;.cr 

" W J2 ( C Rfe', CT c Lfe, CT + /I.C.), (4) 

fc,fc',<7 

where /3 = L, R. Hl (Wr) describes the left (right) 
superconductor with the s-wave order parameter Ae lSL 
(Ae l6,R ). cpk,cr is the annihilation operator for electrons 
with spin a and wave vector k in the superconductor 
(3. Hb represents the QD with the discrete energy level 
£d and the Coulomb interaction U, where d a annihilates 
an electron in the QD and n a — S a d a . The electron 
tunneling is described by Ht, in which t is the tunneling 
amplitude between the superconductors and the QD, and 
w is the direct tunneling amplitude across the insulating 
barrier. 

The Josephson current is simply given by 



e{N L 



It + Iw , 



(5) 



(t,t')= (9) 



(cU A (t')d\(t)) (cp^(t')d\{t)) 



In equilibrium state, the lesser Green function is ex- 
pressed with the advanced and retarded Green functions 
as 



>) = /H (g?>) - g[>)) , (io) 



where f(u>) is the Fermi distribution function. Gfj 

(G^j) is the advanced (retarded) Green function and 

i,j = {Lfc,Rfc,d}. Since (G^.)t = G r ]A , all we have to 
calculate is the retarded Green function. With use of the 
spinor field operators 
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we express the retarded Green function as 



G[>) = ((V^j)L- 



(11) 



(12) 



In calculating the Green functions, we employ the 
equation of motion method. The equation of motion is 
given by 

w«^j»o, = ({«]}> + «^; [H,$])) u . (13) 

It is straightforward to find that all the Green functions 
are expressed in terms of the dot Green function as 
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where we have taken the summation over wave vectors in 
the right hand sides. Here L = R, R = L and T — —ta z . 
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Gf ] r and G^ r are the Green functions for T = which 
are expressed as 



G 



(0)r 



G /3/3 



c^r 1 - 



(17) 
(18) 



where W = —wa z . The unperturbed Green function of 
the superconductor /3 is given by 



VA 2 



-uj Ae 1 ^ 
" V Ae _l9 ' 3 -w 



(19) 



where ^ is the normal density of states at the Fermi level. 
Then, what we have to know in calculating the Josephson 
current is only the dot Green function G dd . 

The Dyson equation for the dot Green function is writ- 
ten as 



G r dd = g r d + G r dd V r T g r d + D r Ug r d , 
where U = Ua z , and 
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(20) 



(21) 



(22) 



Here D t (<jj) — ((V'dj V d ))u> is the higher order Green func- 
tion with 



n T 4 



In order to solve Eq. (|20|) . we decouple D as 



(23) 



((^^))^((^^( ( ^} ) " ( ( ^ T> ), (24) 

where the matrix elements are given by solving the usual 
self-consistent equations 



1 

2tt7 



( 6 U")) n , (25) 

(rf M> = ^/M^K) 21 ' (26) 

and the other elements obey similar equations. The de- 
coupling Eq. (|24[) corresponds to the usual HFA. Sub- 
stituting Eq. (|24|) into (|20|) . the dot Green function is 
obtained as 



Gr 



(1<1 



(27) 



where g~ is the dressed Green function which is expressed 
as 

{9 *> ~ { -U(d\d\) u; + e d + U{n r ) J' (28) 
Then, we can finally estimate the Josephson current. 




FIG. 2: The critical currents in the non-interacting case 
(U = 0) are plotted against the QD energy level £d- The 
transmission probability r = 4£/(l + £) 2 is changed from 
(bottom) to 1 (top) with a step 0.1. The coupling strength is 
chosen to be T/A = 0.1, which is adapted to all figures below. 
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FIG. 3: (a) ABSs below the Fermi level for r = 0.3. The full 
(dotted) line corresponds to the primary (secondary) ABS, 
i.e., lua = wi (1^2). (b) Current contributions from each of 
the ABSs and the continuous spectrum for r = 0.3. The 
dashed line indicates the continuous spectrum contribution. 



III. RESULTS 

The Josephson current is obtained as a function of the 
phase difference 8 = #l — 9r and the QD energy level 
e d . We define the coupling strength between the QD and 
each of the superconductors by T = nut 2 . The direct 
tunneling across the insulating barrier is characterized 
by £ = ir 2 v 2 w 2 . Because the current-phase relation is 
roughly given by J = I c sin#, for simplicity, we determine 
the Josephson critical current by I c = 1(9 — tt/2) in this 
paper. All the calculations below are done assuming that 
the system is at the absolute zero of temperature. We set 
T/A — 0.1 throughout this paper. 
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The Josephson current consists of the two contribu- 
tions; one from the discrete Andreev bound states (ABSs) 
inside the gap \ui\ < A and the other from the continuous 
spectrum outside the gap \cu\ > A. The ABSs oja are de- 
termined as the poles of the dot Green function, whereas 
the poles of the other Green functions such as G^ r do 
not contribute to the current. 



A. Non-interacting case (U = 0) 

We show the dependence of the critical current on the 
QD energy level in Fig. [2] These line shapes remind us 
of the Fano line shapes which are seen in the differential 
conductance of normal metal systems. Far from the res- 
onance (|ed| — > oo), the currents asymptotically go to a 
certain value, which is given by 



h = 



-:A 



Tsmf 



(29) 



2^ 1 



rsm 



where r = 4£/(l + £) 2 is the transmission probability 
through the direct channel. 

The ABSs below the Fermi level are shown in Fig.[3fa). 
One can see that there are the two ABSs, which we call 
"primary ABS" uj\ and "secondary ABS" loi respectively. 
The secondary ABS is present in the region 



r a r 

^ - A < e d < -= 



(30) 



and disappears for £ = 0, i.e., in the S/QD/S junction. 
We see that the ABSs represent the assymetric depen- 
dence on the QD energy level being affected by the inter- 
ference. FigureJ^b) shows the current contributions from 
each of the ABSs and the continuous spectrum outside 
the gap. It should be noticed that the contribution from 
the primary ABS becomes negative at a certain region, 
whereas the sum of the three is always positive. 

In the previous work, it was claimed that 0-ir transition 
occurred around the dip even in the case U = 0^2. We 
however cannot see the transition in the whole parame- 
ter regime. This difference is caused by the fact that the 
contribution from the continuous spectrum was not taken 
into account. It is known that the contribution from 
the continuous spectrum is not zero in S/QD/S junc- 
tions^^. If we choose the appropriate parameter values 
( r/A = 4, t = 0, 0.6, 1 ), the contribution from the pri- 
mary ABS shows almost the same behavior as Fig. 5(a) 
in Ref. []J. 



B. Interacting case (U 7^ 0) 

We show in Fig. [4] the critical currents as a function of 
the QD energy level at U/A — 10 for different values of 
r. The full line in Fig.fJJa), which is for r = 0, represents 
the current flowing through the S/QD/S junction. As is 
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FIG. 4: Critical currents as a function of the QD energy level 
at U/A = 10 for (a) r = (full line), r = 0.0005 (dashed line), 
r = 0.0015 (dotted line) (b) r = 0.002 (full line), r = 0.003 
(dashed line), r = 0.005 (dotted line) and (c) r = 0.06 (full 
line), r = 0.08 (dashed line), r = 0.1 (dotted line). The inset 
of (b) is an extended figure around the resonance = — U. 



verified by the previous theoretica l 27 ! 28 and experimen- 
tal^— work, a negative critical current (i.e., 7r-junction) 
is found in the magnetic region — U < < 0, in which 
the average electron number in the QD is close to unity. 

With increasing r, the peak structure gradually 
changes because of the interference effect (Fano effect). 
In Figs. QJa) and (b), focusing on the region < —U, 
the critical current near the resonance = ~U decreases 
as t increases. Then, it becomes negative at a certain 
value of r, which means that 0-7T transition occurs. We 
can clearly see in Fig. [5] that the current-phase relation 
changes from / = l/clsinf? to / = — | J c | sin 6 1 by sweep- 
ing ed across the transition point. Note that this 0-7T 
transition occurs in the non-magnetic region < — Z7, 
where the average electron number is close to two. Thus, 
the transition is induced not by the same mechanism as 
that in S/QD/S junctions. We find in Fig-HJc) that the 
critical current is always positive for large r. 

Figure [6] shows the critical currents against the QD 
energy level for different interaction strengths U/A. For 
small U/A, the peak structure around the resonance £d = 
— U drastically changes with U/A. Meanwhile, for large 
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FIG. 5: Current-phase relation at U/ A = 10 for different e d 
which are marked by the filled circles in the inset of Fig.[4](b). 
Explicitly e d /A = -14.5 (full line), e d /A = -13.5 (dashed 
line), e d /A = -12.1 (dotted line), e d /A = -11 (dot-dashed 
line) . 



U/A, the peak structure does not depend on the value 
of U/A (it just shifts the position of the resonance.). 
Here the negative critical current in the dip region of the 
resonance e d = — U is found again for large U/A. 

As is shown above, the Fano effect induces 0-ir transi- 
tion with large U /A and appropriate values of r. Almost 
all the behaviors shown in this section are also seen in 
the results obtained by the Hubbard-I approximation 29 , 
which we show in the Appendix. 
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FIG. 6: Critical currents as a function of the QD energy level 
for r = 0.02 with (a) U/A = 1 (full line), U/A = 1.25 (dashed 
line), U/A = 1.5 (dotted line), U/A = 1.75 (dot-dashed line) 
(b) U/A = 2 (full line), U/A = 3 (dashed line), U/A = 4 
(dotted line), U/A = 5 (dot-dashed line) and (c) U/A = 6 
(full line), U/A = 8 (dashed line), U/A = 10 (dotted line), 
U/A = 12 (dot-dashed line) 



IV. PERTURB ATIVE APPROACH 



It is important to find the physical origin of the 0-ir 
transition and the parameter regime where the system 
is a 7r-junction. To do this, we recalculate the Joseph- 
son current by a simple perturbation theory with re- 
spect to the electron tunneling. The Josephson current 
flowing through our system has three components, I = 
h + ^qd + iint ■ Here ij and Jqd are the non-interference 
components associated with the path through the insu- 
lating barrier and the QD respectively, and Ji n t is the in- 
terference component. They are expanded as = l w i + 

i«j 4 H , iQD = if « +i*8 H , and l mt = I wt a +I w 2 t 2 -\ , 

where I w n t m is the term which is proportional to w n t m . 
We calculate the lowest order terms of each of the com- 



ponents. Here we show the results only for e d < — U: 
3 A /w\ 2 x-^ 2A 3 



~7T lAV ^ 

k,k' 



' E k E k >(E k + E k >) 



(31) 



2A 5 



h \AJ ^E k E k ,{E k -e d -U){E k , -e d -U) 



1 



E k +E k , 2e d + U, 

2 2A 4 



(32) 



l Wt 2 ,C 



•t (I) (I) E 



k.k 

2 1 



E k E k ,{E k , -e d -U) 



E k +E k , E k -e d -UJ' 



(33) 



where I w ^ t "\c = Iw n t m \e=-rr/2> and E k = + A 2 . One 
can see that only the term of the interference component 
is negative. It indicates that the interference effect can 
induce 0-7r transition when | I wt 2 c | becomes large enough. 
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In fact, I w tz, c changes its sign at the resonances and the 
particle-hole symmetric point e d = — [7/2; 



<0 (e d <-U) 

> (-U <e d < 

< (-17/2 < e d < 0) 

>0 (0<e d ) 



(34) 



This sign change leads to the characteristic peak-dip 
structure of the Fano resonance. 

The term I wt z c originates from the tunneling processes 
in which the two electrons of a Cooper pair take the two 
different shortest paths. Two of such tunneling processes 
are schematically represented in Fig. [7] In each of the 
tunneling processes, the two electrons acquire the differ- 
ent additional phases 4>qb and (f)\, respectively. Then, 
the phase of the transported Cooper pair is shifted by 
4> = 4>i + </>qd- The negative sign of I w t 2 , c is attributed 
to this phase shift. 

To examine what value <p takes in each of the tunnel- 
ing processes, let us consider two of them. Figure [7^a) 
shows one for the case e d < —U, where the QD energy 
level is almost doubly occupied. In this case, to transfer 
an electron from the left to right superconductor through 
the QD, first an electron in the QD must tunnel out to 
the right superconductor. Then an electron in the left su- 
perconductor tunnels into the QD. Meanwhile in the case 
— £7 < e d < — [7/2 as in Fig.[7Jb), it is energetically favor- 
able that first an electron tunnels into the QD. Compar- 
ing these two tunneling processes, we see that the order of 
the elementary tunneling events is permuted. Because of 
the anticommutation relation of electron operators, this 
permutation results in the 7r-difference of the additional 
phase 0qd- 

Calculating the tunneling amplitudes actually, one 
finds that 0qd = (tt) for the tunneling process rep- 
resented in Fig.[7ja) (Fig.[7[b)). On the other hand, the 
electrons tunneling through the insulating barrier get the 
additional phase fa = Arg(— w) = it. As a result, the 
phase shift of the transported Cooper pair is given by 
4> = 7r (0) for the tunneling process shown in Fig. [7ja) 
(Fig.Etb)). 

The phase shift <j) — 7r makes the tunneling process con- 
tribute to the critical current negatively. For e d < — U, 
<f> = 7r in all the tunneling processes which contribute to 
I w t 2 ■ This is the main origin of the Fano-induced 0-7r 
transition in our system. Simply put, the essential origin 
of this transition is that the two electrons of a Cooper pair 
can take different paths and acquire additional phases the 
sum of which equals ir. If one changes the sign of w, the 
e d dependence of the Josephson current is converted as 
7(e d ) — * J(— e d — U) because of the particle-hole symme- 
try. However, one can always find the negative critical 
current in a non-magnetic region independently of the 
sign of w. 

In the limit e d —>—[/ — + , where \I w t^ c | has a maxi- 
mum value, we can perform the summations in Eqs. (|31|) 







(b) 







FIG. 7: Schematic representation of tunneling processes 
which contribute to I wt 2 for (a) ea < —U, and (b) — U < 
e c \ < —U/2. Dashed line in the energy diagrams of the QD 
means the Fermi level chosen as the zero of energy. 



(|3"3")l analytically; 



J 2 -^2C 
h 

_ eA 2r 2 
It4 < c ~ Itf 



2A 

IT 



*-wt 2 ,c 



(35) 
(36) 

(37) 



Substituting these expressions to I c w l w i c + I t 4 c + 
Iwt 2 .c < 0, we estimate the conditions for the negative 
critical current at e d — * — U — + ; 



C < C < Cc + , 
u > u c , 



(38) 
(39) 
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where 




1 - 2/tt 



2.75A. 



(40) 
(41) 



From Eqs.(|38|) and (j39|) . we see that the negative crit- 
ical current is possible only for large U and the appro- 
priate values of £. It is consistent with the results ob- 
tained by the HFA shown in the previous section. The 
critical values of r for U/A = 10 and T/A = 0.1 are 
given by = 4(±/(l + C^) 2 ~ 0.003, 0.112. These 
are somewhat larger than the values Tjjfa c « 0.002,0.05 
estimated from the results of the HFA. Eq. (|4U|) tells us 
that w/A ~ (t/A) 2 is necessary for the negative crit- 
ical current. In this regime, the interference and non- 
interference components are the same order, and the 
equation / = l w i + I t i + l wt i is correct within the second 
order of T/A. 



the HIA is more reliable at least near the resonances 
where the charge fluctuation plays an important role. 

First we derive the equation of motion for D in 
Eq. (f2"0"|) , and then decouple the higher order Green func- 
tions which arise from ({ipa; [H, <p d ]))u- After that, we 
get the following equation; 



+((< t p d ;4)) u kfg r u , (Al) 



where 
5u : 



Ilj - <r z (e d + U) 

,_+_+. 



* ! ' ±<44) (n T ) ; ' 



V. SUMMARY 

We have investigated the Fano effect in a Joscphson 
junction including a QD in the weak coupling regime 
Tk <C A. To treat the Coulomb interaction inside the 
QD, we employ the HFA. The interference effect between 
the direct and resonant tunneling processes causes the 
characteristic Fano line shape. Furthermore, it is found 
that our system behaves as a 7r-junction even in the non- 
magnetic region ed < — U for large U/A. These results 
suggest that 0-ir transition is induced by the cooperation 
of the Fano effect and the Coulomb interaction inside the 
QD. 

A superconducting quantum interference device 
(SQUID) including carbon nanotube quantum dots was 
already achieved in the recent experiment 2 ^. We expect 
that our theory will be experimentally verified using the 
similar system in the near future. 
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APPENDIX A: HUBBARD-I APPROXIMATION 

In this Appendix, we calculate the Josephson current 
by the Hubbard-I approximation 2 ^ (HIA) generalized to 
the superconducting case. By using the HIA, one can 
take into account the correlation (Sn-[Sni) beyond the 
simple HFA. Since this correlation largely contributes to 
the charge fluctuation-^/ ((Sn^ + 5ni) 2 ), we believe that 



Equation (IA4[) seems complex to treat for a first look. 
However, assuming that the current through the QD is 
uniform, one can see that the following relation holds 



= 0. (A5) 



In addition, we assume that one can disregard the time 
derivative of the pair amplitude inside the QD, i.e., 



(^Md\)) = (2e d + U)(d\d\) 



(It 



-t 



£(<4,X> + <4 c k;>)=°> ( A6 ) 



which guarantees the condition (GJJd)^ 
Eq. (|A4| is simplified as 



G a dA . Thus 



A 



(2e d + U) 





t"l/ 



{d\d} 







(A7) 



The dot Green function is formally obtained by substi- 
tuting Eq. (|AT|) to (J20J) as 



G. 



dd 



(9 r uY 



A+U 

- t 



UK. 



ATU 



(A8) 



From this dot Green function, we estimate the Josephson 
current in the same way as in Sec. II. 

In Fig. [8l we show the results of the HIA together with 
the ones of the HFA. In the case t ^ 0, the two results 
show almost the same behavior, and the negative critical 
current in the region ed < — U is found in both the results. 
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Since, compared to the HIA, the charge fluctuation is un- 
derestimated in the HFA, the difference between the two 
results becomes large as ed approaches the resonances. 

Meanwhile, for r = 0, the obvious difference is seen 
in the region — U < < 0. That is, a negative critical 
current is not seen in the results of the HIA. It means that 
the HIA fails to demonstrate the 7r-j unction behavior in 
the S/QD/S junction. This is because, in the HIA, one 
overestimates the spin fluctuation which is considered to 
be small in this region. Hence, the spin doublet state is 
not well described. This discrepancy should be corrected, 
but this is out of the scope of this paper. 



FIG. 8: Plots of the critical currents versus the QD energy 
level at [//A = 10 for (a) r = 0.01 and (b) r = 0. The 
full (dotted) lines represent the results obtained by the HFA 
(HIA). 



* Electronic address: osawa@kh.phys.waseda.ac.jp 

1 T. Hayashi, T. Fujisawa, H. D. Cheong, Y. H. Jeong, and 
Y. Hirayama, Phys. Rev. Lett. 91, 226804 (2003). 

2 D. L. Huffaker, G. Park, Z. Zou, O. B. Shchekin, and D. 
G. Deppe, Appl. Phys. Lett. 73, 2564 (1998). 

3 L. Zhuang, L. Guo, and S. Y. Chou, Appl. Phys. Lett. 72, 
1205 (1998). 

4 J. Kondo, Prog. Theor. Phys. 32, 37 (1964). 

5 S. M. Cronenwett, T. H. Oosterkamp, and L. P. Kouwen- 
hoven, Science 281, 540 (1998). 

6 F. Simmel, R. H. Blick, J. P. Kotthaus, W. Wegscheider, 
and M. Bichler, Phys. Rev. Lett. 83, 804 (1999). 

7 M. Eto and Y. V. Nazarov, Phys. Rev. Lett. 85, 1306 
(2000). 

8 U. Fano, Phys. Rev. 124, 1866 (1961). 

9 R. Schuster, E. Buks, M. Heiblum, D. Mahalu, V. Uman- 
sky, and H. Shtrikman Nature 385, 417 (1997). 

10 H. Lu, R. Lii, and B. F. Zhu, Phys. Rev. B 71, 235320 
(2005). 

11 A. Ueda and M. Eto, Phys. Rev. B 73, 235353 (2006). 

12 K. Kobayashi, H. Aikawa, S. Katsumoto, and Y. Iye, Phys. 
Rev. Lett. 88, 256806 (2002). 

13 Z.Y Zhang, J. Phys. :Condens. Matter 17 4637 (2005). 

14 S.-G. C heng, Y. Xing, X. C. Xie, and Q.-F. Sun, 
larXiv:0806.2176l 

15 E. Vecino, A. Martm-Rodero, and A. L. Yeyati, Phys. Rev. 
B 68, 035105 (2003). 



16 M. S. Choi , M. Lee, K. Rang, and W. Belzig, Phys. Rev. 
B 70, 020502(R) (2004). 

17 F. Siano and R. Egger, Phys. Rev. Lett. 93, 047002 (2004). 

18 K. Grve-Rasmussen, H. I. j0rgensen, and P. E. Lindelof, 
New J. Phys. 9, 124 (2007). 

19 C. Buizert, A. Oiwa, K. Shibata, K. Hirakawa, and S. 
Tarucha, Phys. Rev. Lett. 99, 136806 (2007). 

20 C. Karrasch, A. Oguri, and V. Meden, Phys. Rev. B 77, 
024517 (2008). 

21 J.-P. Cleuziou, W. Wernsdorfer, V. Bouchiat, T. Ondar- 
cuhu, and M. Monthioux Nature Nanotech. 1, 53 (2006). 

22 R. Lopez , M.-S. Choi, and R. Aguado, Phys. Rev. B 75, 
045132 (2007). 

23 H. Pan, and T.-H. Lin, J. Phys. :Condens. Matter 17 5207 
(2005). 

24 G. Sellier, T. Kopp, J. Kroha, and Y. S. Barash Phys. Rev. 
B 72, 174502 (2005). 

25 J. A. van Dam, Yuli V. Nazarov, Erik P. A. M. Bakkers, 
Silvano De Franceschi, Leo P. Kouwenhoven, Nature. 442, 
667 (2006). 

26 A. A. Clerk and V. Ambegaokar, Phys. Rev. B 61, 9109 
(2000). 

27 D. Matsumoto, J. Phys. Soc. Jpn. 70, 492 (2001). 

28 Y. Zhu, O. F. Sun, and T. H. Lin, J. Phys.:Condens. Mat- 
ter, 13 8783 (2001). 

29 C. M. Varma and Y. Yafet Phys. Rev. B 13, 2950 (1976). 



